Modelling COVID-19 pandemic control strategies in metropolitan and rural health districts in New South Wales, Australia

COVID-19 remains a significant public health problem in New South Wales, Australia. Although the NSW government is employing various control policies, more specific and compelling interventions are needed to control the spread of COVID-19. This paper presents a modified SEIR-X model based on a nonlinear ordinary differential equations system that considers the transmission routes from asymptomatic (Exposed) and symptomatic (Mild and Critical) individuals. The model is fitted to the corresponding cumulative number of cases in metropolitan and rural health districts of NSW reported by the Health Department and parameterised using the least-squares method. The basic reproduction number \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$({\mathrm{R}}_{0})$$\end{document}(R0), which measures the possible spread of COVID-19 in a population, is computed using the next generation operator method. Sensitivity analysis of the model parameters reveals that the transmission rate had an enormous influence on \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathrm{R}}_{0}$$\end{document}R0, which may be an option for controlling this disease. Two time-dependent control strategies, namely preventive (it refers to effort at inhibiting the virus transmission and prevention of case development from Exposed, Mild, Critical, Non-hospitalised and Hospitalised population) and management (it refers to enhance the management of Non-hospitalised and Hospitalised individuals who are infected by COVID-19) measures, are considered to mitigate this disease’s dynamics using Pontryagin’s maximum principle. The most sensible control strategy is determined through the cost-effectiveness analysis for the metropolitan and rural health districts of NSW. Our findings suggest that of the single intervention strategies, enhanced preventive strategy is more cost-effective than management control strategy, as it promptly reduces COVID-19 cases in NSW. In addition, combining preventive and management interventions simultaneously is found to be the most cost-effective. Alternative policies can be implemented to control COVID-19 depending on the policymakers’ decisions. Numerical simulations of the overall system are performed to demonstrate the theoretical outcomes.

It can be observed that the cumulative cases and incidences rates are higher in metropolitan health districts compared to rural districts.
Concerns about COVID-19 occurring in places of positive test results have been made real from numerous outbreaks that have already taken place throughout the pandemic 3 . Many countries have responded to the threat by diverting persons away from COVID-19 test centres. However, peer-reviewed literature regarding COVID- 19 and places of positive test results is restricted to opinion pieces or commentaries. These universally emphasize the vulnerabilities of COVID-19 test centres during this pandemic. However, they make no attempt to enumerate the potential cause and consequences of an outbreak, nor set targets for inmate reductions. COVID-19 outbreaks in COVID-19 test centres (e.g. in China, Canada, and the US) have largely gone unremarked in the peer-reviewed literature. In contrast, the 'super spreader' event onboard the Diamond Princess cruise ship has been the subject of numerous epidemiological publications 4 leading to many lessons learnt.
Experience with respiratory outbreaks and other diseases have repeatedly demonstrated that no matter how isolated a community is, it is not necessarily insulated from infection. Even where there are no complicating factors, such as the age distribution or the presence of individuals with greater susceptibility in the enclosed population, their organization tends to increase transmission and secondary infection risk 5 . In this way, conditions of positive test results have been described as 'amplifiers of infectious disease' 6 . Numerous outbreaks, particularly reports of influenza 7 and tuberculosis 8 demonstrate the propensity for outbreaks to spread rapidly.
Modelling has been utilised as a tool to address gaps in knowledge and inform detention health policies in the prevention and control of infectious disease in detention settings 9 . Researchers have been using various models to provide insights and guide interventions to control the spread of COVID- 19. These have been undertaken in various settings such as the epicentre of Wuhan, China 10 , in high-income countries to predict health service needs 11,12 , in middle-income countries to check the effects of total population and population densities 13,14 to refugee camps and other low-income settings 15,16 .
On the 10 th of February 2022, 988,357 out of the 2,433,278 COVID-19 cumulative cases (around 40%) in Australia occurred in NSW 17 . Strong contingency planning is needed so that health authorities can act accordingly 18 . To this end, mathematical modelling offers an alternative source of information to enlighten preparedness and preventative interventions. Given the differences between metropolitan and rural areas in NSW, a model that allows different parameters under different control settings and locations is needed. Motivated by the continuing and distinctive spread of COVID-19 infections in metropolitan and rural areas of NSW (see Fig. 2), this study is aimed at developing a model that can be used to represent the outbreak of COVID-19 in NSW, hence providing practical recommendations for health and justice authorities to support preparedness measures to better protect the vulnerable population and surrounding communities.
We develop a novel nonlinear deterministic SEIR-X model, or the SEMCNHRD model, for COVID-19 incorporating transmission routes from two infectious classes, including mild and critical individuals, which is parameterized and analyzed based on the cumulative number of reported cases. The model behaviour over a long period is examined qualitatively and quantitatively. Sensitivity analysis is used to assess the impact of fluctuations in model parameters for the epidemiological threshold on the formulation of strategies required to control the spread of the disease. Further, the model can be used to derive the optimal control strategy and the optimal levels of time-dependent preventive and management procedures to be implemented in metropolitan and rural districts to diminish the number of cases in the population effectively. The cost-effective intervention, which combines both preventive and management control strategies, capable of flattening the number of cases curve over a fixed time interval, is recommended and compared between metropolitan and rural health districts of NSW.
The remainder of this paper is structured as follows. In "Methods and materials" section, we present the SEIR-X model with differential infectivity, the model fitting procedures and the sensitivity analysis of the model outputs. "Optimal control strategy and analysis" section reports the optimal control strategy and analysis results. Finally, in "Discussion", we provide a summary of significant findings, discuss their importance for public health policy making and propose guidelines for future efforts.

Methods and materials
Model development. The following model, which we call the SEIR-X model, is proposed for modelling the transmission of COVID-19. The model consists of the following mutually exclusive compartments: Susceptible S(t) , uninfected individuals who are susceptible to the COVID-19 infection; Exposed E(t) , those who are infected but have not yet entered the active COVID-19 stage; M(t) , Mild individuals who are infected, infectious and have mild respiratory illness symptoms such as nasal congestion, runny nose and a sore throat; C(t) , Critical individuals who are infected, infectious and have severe symptoms including shortness of breath, chest discomfort and bluish face. After seeking medical advice, critical individuals are classified as either N H -nonhospitalised individuals but still infected, and H(t)-hospitalised individuals who are still infected. The last two compartments are the Recovered R(t) , who were previously infected and were successfully treated, and Death D(t) . At the same time t, an individual is classified into one and only one compartment. The total population size N(t) is assumed to be a constant at time t and well mixed: To ensure the population size remains constant, we replace all deaths by newborns in the susceptible class. This includes deaths through natural causes, which occur in all states at a constant rate µ , and COVID-19 related deaths which occur at a constant rate ω . Susceptible individuals may be infected with a circulating strain of COVID-19 at the rate = β(M(t) + C(t)) and move to the corresponding exposed class E(t) . Here, β is the probability of susceptible individual contracts infection after contact with Mild or Critical individuals with COVID-19. Those with latent infection progress to mild (the M(t) class) due to reactivation of the latent www.nature.com/scientificreports/ infection at an average period α . However, some Mild individuals move to the recovery class R(t) at an average period ρ due to natural recoveries and the rest of the Mild class individuals move to the critical compartment at an average period φ due to the progression and possibly comorbidities with other diseases, including hypertension, diabetes, cardiovascular disease, and respiratory system disease 19 . A proportion (type of ratio) of the Critical individuals move to the Non-hospitalised and Hospital compartments at an average period γ 1 and γ 2 , respectively. Some Non-hospitalised individuals progress to the recovered compartment R(t) at an average period τ 1 through treatment and the rest progress to the death compartment (D) at an average period δ 1 . Similarly, some of the Hospitalised individuals progress to the recovered compartment R(t) at an average period τ 2 through treatment and the rest progress to the death compartment (D) at an average period δ 2 . A flow diagram of our proposed model is presented in Fig. 3.
In this case, the model can be expressed by the following deterministic system of nonlinear ordinary differential equations: The SEIR-X (SEMCNHRD) model structure: the population is divided into the following eight classes: susceptible, exposed (and not yet symptomatic), infectious (symptomatic) i.e. mild (mild or moderate symptom) and critical (severe symptom), death and recovered (i.e. isolation, recovered, or otherwise noninfectious). www.nature.com/scientificreports/ Given the non-negative initial conditions for the system above, it is straightforward to show that each of the state variables remains non-negative for all t > 0. Moreover, summing Eqs. (2)-(9) we find that the size of the total population, N(t) satisfies i.e. N(t) = constant.
This shows that the total population size N(t) is a constant and it naturally follows that each of the compartment states S, E, M, C, N H , H, RandD is bounded. Given the positivity and boundedness of the system solutions, the feasible region for Eqs. (2)-(9) is given by where D 1 is positively invariant.
Basic reproduction number. The basic reproduction number ( R 0 ) is defined as the expected number of secondary cases created by a single infectious case introduced into a totally susceptible population. The disease can spread in a population only if the basic reproduction number is greater than one. An epidemic occurs when an infection spreads through and infects a significant proportion of a population. A disease-free population is possible when the basic reproduction number is less than one, which means that the disease naturally fadesout 20,21 .
There are eight states in the modelling system in which five belong to the infected states, i.e. E, M, C, N H and H , and three are uninfected states, i.e. S, R and D . At the infection-free steady-state (3)-(7) are closed, in that they do not involve the derivation of S from steady state value. Also, R and D do not appear in Eqs. (3)-(7), and for (E, M, C, N H , H) we have the following equations: Here, these Ordinary Differential Equations (ODEs) in (10)- (14) are referred to as the infection subsystem, as they only describe the production of newly infected individuals and changes in the states of already infected individuals.
By setting x = (E, M, C, N H , H) ′ , where the prime denotes transpose, the infection subsystem can be written in the following form: The matrix T corresponds to the transmission, and the matrix to transitions. All epidemiological events that lead to new infections are incorporated in the model via T and other events via . If the infected states are indicated with i and j with i, j ∈ {1, 2, 3, 4, 5} , then the entry T ij is the rate at which individuals in infected state j give rise to individuals in infected state i . The matrices T and admit the form The next-generation matrix, K , is given by 22 (note the essential minus sign) The dominant eigenvalue of K is the basic reproduction number for COVID-19, which represents the average number of new infections produced by one infected individual. Here, the basic reproduction number can be expressed as: The first term represents the probability of becoming infectious once infected with mean infectious period 1 (φ+ρ+µ) . The second term indicates the expected number of infected individuals due to the comorbidities with other diseases with mean infectious period 1 (φ+ρ+µ)(γ 1 +γ 2 +ω+µ) . We provide detailed analysis, including the existence and stability of the equilibrium points, for the proposed COVID-19 model (2) 9) to the actual reported cases in metropolitan and rural health districts in NSW 23 . In order to parameterize the model, we obtained some of the initial parameter values from literature (see Table 1), and others were estimated from data fitting. The estimation of parameters was carried out using the least squares method which minimises summation of the square errors given by Y t, q − X real 2 subject to the COVID-19 model (2) Figure 4 shows the incidence data of COVID-19 (red dash) and the model fitted curve (blue solid curve).  (9) is to variations in each of its parameters in order to advise intervention strategies that will support in bringing down the infection trajectory. Further, sensitivity analysis will help understand what should be prepared or avoided to mitigate the outbreak of the COVID-19 28,29 . For this purpose, we calculate the partial rank correlation coefficient (PRCCs) 30,31 between each of the model parameters and several output variables (Mild and Critical cases) using a Latin Hypercube Sampling. Specially, a uniform distribution is allocated from half to fourfold baseline value (see Table 1) for each model parameters and assigned 100,000 simulations for each. Figures 5 and 6 display the correlations between Mild and Critical cases of metropolitan and rural health areas and the corresponding parameters β, α, φ, ρ, γ 1 , γ 2 and ω . From Figs. 5, 6 and 7, it is observed that Mild and Critical cases have a strong positive correlation with parameters β (transmission rate) and α (progression rate from E to M) in both metropolitan and rural health areas, implying that increasing β and α will rise Mild and Critical cases. On the other hand, parameters φ, ρ, γ 1 , γ 2 and ω have a negative correlation with Mild and Critical cases, implying that increasing those parameter values will reduce Mild and Critical cases. As discussed in earlier sections, the scale and severity of COVID-19 transmission are directly related to the basic reproduction number R 0 . Here, we assessed the sensitivity indices of the reproduction number R 0 . The indices identify how significant each parameter is to R 0 and thus the COVID-19 transmission dynamics, and recognize which area should be focused in terms of intervention strategies.
Thus, from the explicit formula for R 0 , the analytical expression for the sensitivity indices can be derived as a comparative variation in R 0 when each parameter changes using the following equation 32 :   www.nature.com/scientificreports/ where ϒ R 0 q is the sensitivity index of a differentiable R 0 for any parameter, q. The sensitivity indices for the model (2)-(9) are graphically presented in Fig. 7. It can be perceived that of all the positives indices, the effective contact rate,β , is the highest in both metropolitan and rural areas, and therefore the most sensitive parameter. The value of the sensitive index suggests that an increase (or a decrease) in the value of β will increases (or decrease) R 0 by 100%. However, of all the negative indices displayed in Fig. 7, the recovery rate for the Mild class, denoted by ρ , is the most sensitive parameter in both metropolitan and rural areas. An increase (or a decrease) of the value of ρ will decrease (or increase) R 0 by 85%. Moreover, we observed that the progression rate α is comparatively more sensitive in the metropolitan area than in the rural area.
Additionally, contour plots of R 0 as a function of other parameters are displayed in Fig. 8 to determine how variations in these parameters affect the basic reproduction number, R 0 . Figure 8A shows a decrease in R 0 with increasing recovery rates for both progression rates from Critical (C) class to Non-hospitalised and Hospital classes. Further, Fig. 8B shows a decreasing trend of the basic reproduction number with both transmission rate, β , and recovery rate from M to R, ρ . It is conjectured that R 0 can be brought lower than the threshold of one if efforts are geared towards dropping the contact rate while concurrently improving control and treatment of COVID-19 cases.

Ethical approval. This study is based on aggregated COVID-19 surveillance data from the New South
Wales (NSW) in Australia provided by the NSW government. No confidential information was included because mathematical analyses were performed at the aggregate level. We compiled data from the publicly available website https:// data. nsw. gov. au/ search/ datas et/ ds-nsw-ckan-aefcd e60-3b0c-4bc0-9af1-6fe65 2944e c2/ detai ls?q = .
(a) Metropolitan health district (b) Rural health district

Optimal control strategy and analysis
With consideration of the sensitivity result, control measures that were implemented to other diseases models [31][32][33][34][35][36] were considered to see their effectiveness in controlling the spread of COVID-19 in the population of the metropolitan and rural area. This is executed by introducing two time-dependent control variables u 1 (t) and u 2 (t) , defined as follows: i. u 1 (t) denotes the preventive strategy that is the effort at inhibiting the virus transmission and prevention of case development from Exposed, Mild, Critical, Non-hospitalised and Hospitalised population. This can be reached through public health advocacy for social distancing, good personal hygiene, diagnosis campaigns, wearing face masks in public places, education programs for public health, effective treatment with completion, and protective gear for healthcare workers. Noting that u 1 (t) = 1 indicates the policy effectively protects against infection, while u 1 (t) = 0 denotes the absence of the strategy. ii. u 2 (t) indicates control variable to enhance the management of Non-hospitalised and Hospitalised individuals with a view to ensure the rapid provision of additional treatment or oxygen or mechanical ventilation for Non-hospitalised and Hospitalised individuals with mild and severe COVID-19 symptoms. If u 2 (t) = 1 , then the control strategy is effectively managing the disease, while u 2 (t) = 0 means the absence of the control strategy.
The forward-backward sweep method 37 is used to solve the subsequent optimal control problems 38 . The incremental cost-effectiveness ratio (ICER) is used to determine the intervention strategy for the cost-effective analysis that is the best value for money. We provide detailed analysis of the optimal control of our proposed model in the supplementary materials (see the optimal control analysis section).

Simulation of optimal control and cost-effective analysis.
We implemented the Runge-Kutta fourth order forward and backward method using MATLAB programming language to solve the subsequent optimality system which consists of (Supplementary Eq. 12) and (Supplementary Eq. 16) with the characterization (Supplementary Eq. 18) within the period of [0, 100] days. The weight constants adopted for balancing the objective function (Supplementary Eq. 13) are selected to ensure that no term dominates the other. Therefore, we used equal weight constant for minimising the infectious classes, so that a 1 = a 2 = a 3 = a 4 = a 5 = 1. Under other conditions, the weight constants for determining efforts or cost essential to implement the controls are comparatively different, and outcomes in values for a 6 = 50 and a 7 = 100 , which are consistent with previous research 39 . Details of the numerical procedure for simulating the obtained optimality system are contained 37 . Figure 9 establishes how single preventive measure, u 1 (t) , affects the spread of the COVID-19 in the metropolitan and rural areas in NSW. As shown in Fig. 9a,b, to minimise the objective function (S13), the optimal control u 1 (t) is continued at the maximum level (i.e. 100%) for about 40 days for metropolitan and 25 days for rural health districts before relaxing to the minimum in the final time. Also, as expected, the number of COVID-19 infectious individuals are reduced when control is in place. Figure 10 displays the effects of a single management intervention strategy u 2 (t) on the dynamics of COVID-19 virus infection in the metropolitan and rural health districts. We observed that the u 2 (t) intervention has small impact on the Exposed, Mild and Critical class population while it has significant impact on Non-hospitalised and Hospitalised class population in both areas. Our finding is consistent with the reality because if we only take care of the Non-hospitalised and Hospitalised population, then the Critical class population moves to the Non-hospitalised and Hospitalised class population, and that is why implementing u 2 (t) control strategy has significant impact on Non-hospitalised and Hospitalised class population.
(a) RelaƟonship between R 0 and both γ 1 and γ 2 (b) RelaƟonship between R 0 and both β and ρ   www.nature.com/scientificreports/ Figure 11 shows the implication of combining the two optimal controls in bringing down the total number of infectious people to zero in both metropolitan and rural health areas. It is observed that optimal solution is achieved when preventive control strategy (u 1 ) is strictly adhered to at the maximum level of 100% for 45 days and 30 days, while the management control strategy (u 2 ) of the Non-hospitalised and Hospitalised individuals is at a maximum level of around 25% in metropolitan and rural health areas in NSW. It can be seen that the combination of the two control strategies is significantly more effective in decreasing the spread of the COVID-19 virus compared to the implementation of each control strategy individually. This is consistent with previous modelling studies 30,36,40 .  www.nature.com/scientificreports/ It is essential to determine the most cost-effective strategy among the single and combined control strategies to optimally mitigate the spread of COVID-19 at the possible minimum cost. This is performed by associating the differences among each intervention's costs and outcomes; obtained by estimating the incremental cost-effective  www.nature.com/scientificreports/ ratio (ICER), which is defined as the extra cost per additional intervention outcome. Incrementally, when analysing two or more competing intervention policies, one intervention is associated with the next less effective option. The ICER numerator is given by the total difference in intervention costs, active COVID-19 cases averted costs and averted productivity losses if applicable, between each scenario and baseline. The ICER denominator is the total number of active COVID-19 cases averted. Hence, the following formula obtains the ICER:   . Control profile for preventive ( u 1 (t)) and management ( u 2 (t)) strategies and its effects on the COVID-19 dynamics in metropolitan health districts (left hand side) and rural health districts (right hand side). www.nature.com/scientificreports/ The total cost for each of the single implementation and mutual effort of the optimal control strategy is obtainable from the objective function (S13). The cases averted is invaded by computing the difference between infectious individuals with and without control strategy. Let, S 1 , S 2 and S 12 respectively represent single preventive intervention strategy u 1 (t) , single management control strategy u 2 (t) and the combined effort of the two strategies. Table 2 summarises the ICER for each and the combination of the control variables u 1 (t) and u 2 (t) in increasing order of the total infection averted.
The ICER results for S 1 , S 2 and S 12 are calculated using (16) and shown in Table 2 follows. www.nature.com/scientificreports/ Comparing S 1 , S 2 and S 12 , it is seen that combined intervention S 12 is the most cost-effective which reduces a significant number of COVID-19 cases in both metropolitan and rural health areas compared to S 1 and S 2 individually, while S 2 is the least effective intervention strategy among them.

Discussion
COVID-19 is one of the most pressing public health problems in NSW. Overall, the transmission dynamics and epidemiology of COVID-19 in NSW is not entirely understood. The NSW government initiated various intervention programs to eliminate COVID-19 last year. Although COVID-19 control in NSW has remarkably progressed-accurate contact tracing, availability of free diagnostic and treatment services, the participation of many partners such as community health care agencies and general practitioners, newer diagnostic services, sufficient human resources, sufficient capacity e.g. hospital admission capacity, and guidelines-more effort is required. To reduce COVID-19 incidence and prevent deaths from COVID-19 in metropolitan and rural health districts of NSW, we need to identify the critical factors for developing COVID-19 disease, improve preventive and management strategies, treatment effectiveness, and reduce failure of treatment in infectious individuals.
In this study, we presented a mathematical analysis of transmission dynamics of the COVID-19 outbreak to deliver further understandings into the disease transmission and explore potential prevention and control strategies capable of reducing the disease spread at the metropolitan and rural districts for better health management in NSW. A compartmental mathematical model was formulated by subdividing the host population into Susceptible, Exposed, Mild, Critical, Non-hospitalised, Hospital, Recovered, and Death. The nonlinear mathematical model has been established by considering transmission routes from Exposed, Mild and Critical population individuals.
The formulated model was fitted to the reported data at the metropolitan and rural health districts in NSW, Australia. The least-square method was applied to estimates parameters such as transmission rate (β) and progression rate (α) from Mild to Critical. The other parameters values were determined from the country demographic profile and literature review. We obtained the basic reproduction number (R 0 ) , using the next generation matrix. Sensitivity analysis of the model was performed to find the parameters that drive the spread of the virus infection mostly in the population. It was revealed that transmission rate (β) is the most sensitive parameter, which has a positive correlation with R 0 . It meant that decreasing the transmission will reduce the COVID-19 cases for both metropolitan and rural health geographies. Further, the recovery rate (ρ) is the second most sensitive parameter, negatively correlated with R 0 . It refers to increasing this parameter value will decrease the secondary cases of COVID-19.
We implemented an optimal control approach via Pontryagin's Maximal principle 41 and formulated the optimal strategies for controlling the COVID-19 epidemic at metropolitan and rural areas in NSW. Two different control strategies were considered including preventive strategy (u 1 ) and management of non-hospitalised and hospitalised strategy (u 2 ) . Different settings were examined to measure the cost-effectiveness of both preventive and management control strategies. Between the two-single control strategies, the preventive strategy (u 1 ) is better in cost-effectiveness than the management strategy (u 2 ) which reduce a significant number of COVID-19 cases in both metropolitan and rural health districts of NSW in Australia and similar to the previous study 26 . Therefore, our results suggest that the NSW government should improve preventive control interventions when only one control strategy is used. Naturally, this strategy actively decreases and/or stops the contact between susceptible and infectious individuals of COVID-19. However, combined implementation of preventive and management strategies is the most cost-effective measure for reducing the burden of COVID-19 in both metropolitan and rural districts in NSW.
Optimal control strategies have been applied in other endemic settings to minimise COVID-19 cases and intervention implementation costs. Previous studies show that preventive strategy is the best strategy for the single control strategy implementation to decrease COVID-19 cases and intervention costs 39,42,46 , which is similar to our results. However, our finding also suggests that combining control strategy with health management strategy, including enhanced services facilities and health management is the most effective way to decrease the COVID-19 burden in metropolitan and rural districts in NSW, consistent with previous works 39, [41][42][43][44][45] .
Finally, in NSW, infectious disease surveillance is fully recognised, but the risk of bias cannot be precluded. More accurate data should be put in place to address concerns related to COVID-19. Accurate data leads to better estimation of crucial parameters, and this means our proposed intervention to decision support is datadependent. Hence, local, state and federal level policy-makers need to adjust the possibility of under-reporting bias when investigating our outcomes. Therefore, more accurate data could be included in the model to explore the impact of preventive and management interventions on COVID-19 dynamics in NSW. Of the single intervention strategies, enhanced preventive strategy is cost-effective compared to management control strategy and www.nature.com/scientificreports/ prompt reducing COVID-19 cases in NSW. Yet, combining both preventive and management interventions is found to be even more cost-effective. In summary, the paper provides mathematical modelling and optimal control strategy of COVID-19 in metropolitan and rural health districts of NSW for better pandemic management in Australia. We derived the basic reproduction number and found that it's play an important role in the outbreak of COVID-19 in both districts. Sensitivity analysis also performed to identify the most important risk factors and found that transmission rate had largest influence on COVID-19 prevalence. We adopted optimal control analysis via Pontryagin's Maximal Principle and formulated the optimal control strategies for controlling the COVID-19 outbreak in NSW, Australia. However, in our model we considered the total population size is constant and mixes homogeneously. That means that each individual in a compartment has identical disease susceptibility, infectiousness and transmission frequency with others, which neglects social behaviours of children, younger and older people that are known to be important risk factors in the spread of COVID-19. In addition, in our model individual-level COVID-19 natural history is not incorporated and not justified based on any empirical evidence. Therefore, future research could focus on individual-level data and quantifying any data-oriented reporting bias and modelling appropriate priority-based vaccination coverage.
Two different control strategies including preventive strategy and management control strategy were implemented to measure their cost-effectiveness for improved and consistent health service facilities to the communities. Among the two single-controls, the preventive control strategy is the most cost-effective. Therefore, when single control strategy is used, our results suggest that NSW government should improve preventive control strategy, reducing contact between infectious and susceptible individual. Our principal finding is that the combination of preventive and management control strategies such as the rapid provision of additional services means is the most impactful and cost-effectiveness strategy for reducing the spread of COVID-19 in NSW, Australia. The same modelling framework can be utilised for data obtained in other places around the globe.

Data availability
The datasets produced during the study are available from the corresponding author on reasonable request. All data were compiled from the publicly available website https:// data. nsw. gov. au/ search/ datas et/ ds-nsw-ckan-aefcd e60-3b0c-4bc0-9af1-6fe65 2944e c2/ detai ls?q = .